AD-A128  622 


UNCtASSIFIED 


OPTICAL  TRACKING  INVESTIGATIONS  INTO  INAGE  INFORMATION 
AND  ADAPTIVE  GUIDANCE OJ>  FRANK  J  SEILER  RESEARCH  LAB 
UNITED  STATES  AIR  FORCE  ACADEHV  CO  J  D  LEDBETTER 
APR  83  FJSRL-TR-83-8884  F/G  17/8 


1/1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-  1963-A 


FRANK  J.  SEILER  RESEARCH  LABORATORY 


APRIL  1983 


OPTICAL  TRACKING  INVESTIGATIONS  INTO 


IMAGE  INFORMATION  AND  ADAPTIVE  GUIDANCE 


FINAL  REPORT 


CAPTAIN  TAMES  D.  LEDBETTER 


APPROVED  FOR  PUBLIC  RELEASE; 


PROJECT  2305-F2-65 


DISTRIBUTION  UNLIMITED. 


f-V  4  •  ,  VjS 

Vv.-v  v..  ' 


AIR  FORCE  SYSTEMS  COMMAND 


UNITED  STATES  AIR  FORCE 


8  ^ 


T  O  •  -1 


This  document  was  prepared  by  the  Guidance  ami  Lontrol  Division, 
Directorate  of  Aerospace-Mechanics,  Frank  J.  Seiler  Research  Laboratory, 
United  States  Air  Force  Academy,  Colorado  Springs,  CO.  The  research  was 
conducted  under  Project  Work  Unit  Number  2305-F2-65,  Adaptive  Concepts  in 
Electro-Optical  Tracking  Applications.  Captain  James  R.  Ledbetter  was  the 
Project  Scientist  in  charge  of  the  work. 

When  U.S.  Government  drawings,  specifications  or  other  data  are  used 
for  any  purpose  other  than  a  definitely  related  government  procurement 
operation,  the  government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever,  and  the  fact  that  the  government  may  have 
formulated,  furnished  or  in  any  way  supplied  the  said  drawings, 
specifications  or  other  data  is  not  to  be  regarded  by  implication  or 
otherwise,  as  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation  or  conveying  any  rights  or  permission  to  manufacture,  use  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

Inquiries  concerning  the  technical  content  of  this  document  should  be 
addressed  to  the  Frank  J.  Seiler  Research  Laboratory  (AFSC),  FJSRL/NH, 
USAF  Academy,  Colorado  Springs,  CO  80840.  Phone  AC  303-472-3122. 

This  report  has  been  reviewed  by  the  Commander  and  is  releasable  to 
the  National  Technical  Information  Service  (NTIS).  At  NTIS  it  will  be 
available  to  the  general  public,  including  foreign  nations. 

\ 

This  technical  report  has  been  reviewed  and  is  approved  for 
publication. 


—  1 


JAMES  D.  LEDBETTER,  Captain,  USAF  THEODORE  T.  SAITO,  Lt  Col,  USAF 

ject  Scientist  Director,  Aerospace-Mechanics  Sciences 


WILLIAM  D.  SIURU,  JR 


Commander 


Copies  of  this  rejJort  should  not  be  returned  unless  return  is  required 
by  security  considerations,  contractual  obligations,  or  notice  on  a 
specific  document. 


Printed  in  the  United  States  of  America.  Qualified  requestors  may 
obtain  additional  copies  from  the  Defense  Technical  Information  Center. 
[All  other  should  apply  to:  National  Technical  Information  Service 

6285  Port  Royal  Road 
Springfield,  Virginia  22161] 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dara^Fntered) 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

HE  FORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

FJSRL-TR-8 3-0004 

2  GOVT  ACCESSION  NO. 

3  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

Optical  Tracking  Investigations  Into  Image 
Information  and  Adaptive  Guidance 

S  TYPE  OF  REPORT  6  PERIOD  COVERED 

Final  Report 

Sep  79  -  Aug  82 

6  PERFORMING  OTG.  REPORT  NUMBER 

7  AuThORi'j; 

Captain  James  D,  Ledbetter 

8  CONTRACT  OR  GRANT  NUMBERfaJ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Frank  J„  Seiler  Research  Laboratory 

USAF  Academy 

Colorado  Springs,  CO  80840 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  8  WORK  UNIT  NUMBERS 

2305-F2-65 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

April  1983 

13.  NUMBER  OF  PAGES 

74 

14  MONITORING  AGENCY  NAME  »  ADDRESS  (it  different  from  ControlUnt  Office) 

IS.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15*.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  ( of  this  Report) 


Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  30.  II  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  Identify  by  block  number) 

Image  Processing  Reachable-set  Theory 

Texture  Analysis 

Segmentation 

Optical  Tracking 

Adaptive  Guidance 

,20.  ABSTRACT  (Continue  on  rovers*  side  if  necessary  and  Identify  by  block  number) 

A  two  part  study  into  critical  aspects  of  optical  tracking  systems  is  presented. 
The  basic  question  of  what  comprises  information  in  a  textural  image  is 
addressed  in  a  study  to  investigate  the  application  of  an  erosion  process  to  a 
binary  image,.  The  objective  of  the  erosion  is  to  enhance  the  image  properties 
that  match  the  properties  of  the  structural  element,  or  mask,  used  in  the 
erosion  process.  The  resulting  image  can  be  used  for  further  processing,  but 
the  enhanced  properties  would  now  dominate  other  qualities  or  features  in  the 
Image-  -  ‘  ~  (CONTINUED  ON  REVERSE) _ 

DD  I  mnM73  1473  EOITION  OF  I  NOV  65  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fHTien  Data  Entered) 


SECURITY  CLASSIFICATION  of  tu.«-  PAGEflFhan  Data  Enrarad,' 


TABLE  OF  CONTENTS 


Page  No. 

Preface  .  1 

I.  .  2 

1.0  Introduction  and  Background  .  2 

II.  Texture  Analysis  Investigation  .  4 

2.0  Introduction  . .  4 

2.1  Image  Information  . . . . .  5 


2.2  Erosion  Study  .  6 

2.2.1  Conclusions  . 13 

2.3  Statistics  Study  . J-3 

2.3.1  Conclusions  . 15 

III.  Reachable-Set  Guidance  Techt'que  . . 20 

3.0  Introduction  . 20 

3.1  Encounter  Model  . 21 

3.2  Implementation  of  Conventional  Guidance  . 24 

3.2.1  Proportional  . 24 

3.2.2  Pursuit  . 26 

3.3  Advanced  Guidance  . 27 

3.4  Algorithm  Requirements  . 32 

3.5  Computational  Requirements  . 36 

3.5.1  Trigonometric  Look-Up  Table  . . 37 

3.5.2  Multiple  Processors  . 37 

3.5.3  Future  Refinements  . 38 

3.6  Performance  . 43 

3.6.1  Perfect  Measurements  . 45 

3.6.2  Computational  Delay  . 45 

3.6.3  Measurement  Error  Effects  . 47 

Bibliography  . 52 

Appendix  . 54 


LIST  OF  FIGURES 


9 


Figure  No.  Page  No. 

2.1  Seasat  Image  I  .  6 

2.2  Seasat  Image  2  .  6 

2.3  Erosion  Process  .  7 

2.4  Location  of  Edge  (430,  312)  . 8 

2.5  Location  of  Edge  (127,  122)  .  8 

2.6  Seasat  Image  1  -  10%  Threshold  .  9 

2.7  Seasat  Image  1  -  20%  Threshold  .  9 

2.8  Seasat  Image  1  -  30%  Threshold  .  9 

2.9  One  Erosion  -  10%  Threshold  . 9 

2.10  Two  Erosions  -  10%  Threshold  . 9 

2.11  Three  Erosions  -  10%  Threshold  .  9 

2.12  Covariance  For  Edge  (430,312), — Mask . 10 

2.13  Covariance  For  Edge  (127,122), — Mask . 10 

2.14  Covariance  For  Edge  (430,312), — Mask . 12 

2.15  Covariance  For  Edge  (430,312), — Mask . 12 

2.16  Location  of  Edge  (72,372)  . . ..16 

2.17  S08ELSQ  Operation  on  Figure  2.1  . ,.16 

2.18  SEBELSQ  Operation  on  Figure  2.2  . 16 

2.19  32x32  Statistics  for  Edge  (430,312) . 17 

2.20  16x16  Statistics  for  Edge  (430,312) . 17 

2.21  16x16  Statistics  for  Perpendicular  Path  . 18 

2.22  32x32  Statistics  for  Edge  (72,372)  . 18 

2.23  16x16  Statistics  for  Edge  (72,372)  . 19 

2.24  16x16  Statistics  for  Perpendicular  Path  . ,.19 

3.1  Maneuver  Plane  Definition . 22 

3.2  Missile  Acceleration  Vector  . 23 

3.3  Proportional  Guidance  . 25 

3.4  Pursuit  Guidance  . 26 

3.5  Required  Missile  Velocity  Change  . 30 

3.6  Missile  Acceleration  Function  . 31 


ii 


LIST  OF  F IGIJRF.S  (Continued) 


Figure  No.  Page  No. 

3.7  Flowchart  . 33 

3.8  Single  Processor  Solutions  . 40 

3.0  Dual  Processor  Solution  . 41 

3.10  Additional  Uses  For  Second  Processor  . 42 

3.11a  Scenario  2:  Rear  Attack  . 44 

3.11b  Scenario  5:  Side  Attack  . 44 

3.12  Performance  versus  Computation  Delay  . 46 

3.13  Performance  versus  Range  Error  . 48 

3.14  Performance  versus  Velocity  Error . 48 

3.15  Performance  versus  Target  Heading  Error  . 50 

3.16  Performance  versus  Target  Position  Angle  Error  . 50 


PREFACE 


The  work  summarized  in  thi  s  report  falls  into  two  separate,  though 
related,  areas.  Image  processing  techniques  to  extract  target  data  from 
an  optical  sensor  and  the  subsequent  guidance  algorithm  to  use  that  data 
certainly  form  a  closed  loop  system  with  regards  to  purpose.  However,  the 
intent  here  is  to  look  at  basic  concepts  in  both  problem  areas  without 
regard  to  "closing  the  loop"  and  forming  a  complete  tracking  system.  The 
work  reported  on  texture  analvsis  of  image  data  was  done  by  Captain  James 
Ledbetter,  Frank  J.  Seiler  Research  Laboratory.  The  work  on  the 
reachable-set  guidance  evaluation  was  performed  and  written  by  T)r.  Michael 
Larimore  and  Dr.  Claude  Wtatrowski,  University  of  Colorado  at  Colorado 
Springs,  with  minor  revisions  by  Captain  Ledbetter  for  inclusion  in  this 
report. 
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l .0  Introduction  and  Background 

The  continued  advancement  of  solid  state  Imaging  technology  has 

resulted  In  an  Increased  emphasis  on  the  development  of  low-cost,  light 

weight  sensor  arrays  for  Air  Force  surveillance,  reconnaissance,  and 

weapons  control  systems.  The  use  of  these  sensors  as  Imaging  devices  Is 

attractive  due  to  elimination  of  high  voltage  vacuum  tube  technology  of 

1  2 

conventional  vldicon  tubes  ’  .  The  sensor's  small  size  and  low  power 

requirements  Insure  both  linear  and  array  devices  will  have  a  strong 

impact  In  tracking  applications  where  limited  space  and  power  availability 

are  factors.  Additional  technology  verification  has  been  provided  by  an 

Air  Force  Avionics  Laboratory  program  which  conducted  a  parametric 

analysis  of  an  array  tracker  system  to  determine  Its  performance 

characteristics  and  found  that  the  charge  coupled  device  (CCD)  Imaging 

3 

array  was  not  the  limiting  factor  in  most  tracking  applications  .  More 

recently,  the  Jet  Propulsion  Laboratory  has  reported  the  design  of  a  CCD 

4 

tracker  for  a  space  mission  . 

This  report  documents  efforts  undertaken  to  investigate  two  different 
problem  areas  In  optical  tracking.  The  first  area  Is  that  of  target 
segmentation  In  the  optical  field  of  view.  The  segmentation  process  must 
be  performed  early  in  the  processing  of  the  Image  data  and  is  a 
combination  of  boundary  detection  and  texture  analysis  techniques^.  Of 
particular  interest  in  this  study  is  the  information  gained  from  the  image 
texture.  This  interest  is  a  result  of  natural  terrain  where  many  military 
targets  would  be  located.  Such  terrain  is  not  predominately  characterized 
by  tonal  edges  but  by  textural  changes. 

The  second  area  of  interest  is  adaptive  guidance  and  control.  An 
"intelligent"  weapons  system  utilizing  an  optical  sensor  might  ultimately 
incorporate  software  that  would  allow  it  to  adapt  to  the  changing 
scenarios  it  "sees"  by  anticipating  or  predicting  what  a  target  might  do 
to  optimize  its  position  in  an  encounter.  The  present  guidance  techniques 
of  pursuit  or  proportional  navigation  can  be  deceived  by  an  intelligent 
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target  capable  of  radical  manuevers  In  the  final  seconds  of  the 
encounter.  This  report  examines  a  promising  guidance  technique  that  can 
adapt  to  maneuvers.  The  technique  is  evaluated  in  terms  of  its 
implementation  with  present  day  technology. 


SECTION  II.  TEXTURE  ANALYSTS  INVESTIGATION 


2 . 0  Introduction 

Image  segmentation  is  a  maior  component  of  any  machine  image  analysis 
requirement.  Rased  on  the  quality  of  the  segmentation,  other  important 
image  descriptors,  or  features,  can  be  defined  to  further  represent  the 
image  data.  However,  while  humans  find  it  very  easy  to  "see"  a  collection 
of  objects  when  they  view  a  scene,  machines  "see"  only  an  array  of  equally 
weighted  pixel  values  which  vary  in  intensity  based  on  the  amount  of  light 
incident  upon  them.  Obiects,  therefore,  are  not  "seen",  only  pixel 
intensity  from  which  the  objects  comprising  the  image  must  ;  determined. 
Pixel  intensities,  and  the  way  they  are  arranged,  com]  -e  the  hasic 
elements  of  information  available  to  segment  collections  pixels  into 
objects  of  interest  or  regions  which  have  more  or  1  omogeneous 
properties. 

Pixel  intensities  and  their  arrangement  comprise  the  inextricable 
relationship  of  tone  and  texture  properties  in  the  image.  Both  properties 
are  always  present  in  an  image,  one  usually  predominating  over  the  other. 
When  an  area  has  verv  little  variation  in  pixel  intensity,  the  predominant 
property  is  tone.  When  an  area  has  wide  variation  in  pixel  intensity,  the 
predominant  property  is  texture.  The  size  of  the  area  in  this  distinction 
Is  critical.  Its  crucial  nature  arises  when  ^escribing  a  g*ven  texture  in 
terms  of  its  tonal  primitives  and  a  given  spatial  organization.  The  term 
"spatial  organization"  requires  the  declaration  of  the  size  of  the  area 
concerned.  For  that  local  area,  the  texture  Is  then  the  combination  of 
one  or  more  tonal  primitives,  regions  with  tonal  properties,  and  a  spatial 
rule  specifying  their  arrangment.  The  segmentation  of  images  where  the 
information  varies  between  the  tonal  properties  and  the  textural 
properties  is  very  difficult  unless  a  priori  knowledge  is  available  on  the 
statistics  or  properties  of  the  texture.  Without  this  knowledge,  local 
area  operations  must  be  used  to  discern  the  statistics  of  the  texture. 
However,  the  computation  of  these  properties  is  a  function  of  the  size  of 
the  local  area.  This  interrelationship  usually  results  in  the  boundary 
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between  two  textural  areas  being  blurred  or  broadened  due  to  the  averaging 
of  the  properties  of  the  regions  when  the  local  operator  is  straddling  the 
boundary . 

When  regions  of  uniform  texture  are  defined  there  are  many  different 
ways  to  approach  deriving  the  properties  of  the  textures.  A  very  useful 
survey  of  differing  approaches  to  texture  definition  is  given  in  Reference 
6.  One  method  which  has  proven  very  useful  in  classifying  textures  is 
that  of  determining  concurrence  matrices  and  defining  the  texture  based  on 
features  of  those  matrices^.  The  power  of  this  approach  is  that  it 
characterizes  the  spatial  interrelationships  of  the  grey  tones  in  a 
textural  pattern.  Its  weakness  is  that  it  does  not  capture  the  shapes  of 
the  tonal  primitives. 

In  this  section  of  the  report,  information  on  two  investigations  into 
texture  segmentation  techniques  is  presented.  The  image  used  was  almost 
entirely  textural  in  content.  The  primary  Interest  was  in  determining  a 
fast,  simple  "information  indicator"  that  could  be  used  to  enable  an 
adaptive  segmenter  to  switch  between  tonal  and  textural  operations  in 
finding  edges  for  further  boundary  definition.  Section  2.1  describes  the 
image  used.  Sections  2.2  and  2.3  describe  an  erosion  study  and  a 

statistics  study,  respectively,  of  different  textural  edges  in  the  image. 

2.1  Image  Information 

The  images  used  in  these  investigations  were  512x512  subpartitions 

from  a  4000x3000  pixel  image  from  the  Seasat-A  synthetic  aperature  radar 
(SAR)  sensor.  This  sensor  operated  from  July  through  early  October  1978 
and  generated  a  large  amount  of  land  and  sea  data.  Since  it  is  a  radar 
imager,  the  information  In  the  images  i3  almost  entirely  textural  in 

content.  Figures  2.1  and  2,2  are  the  images  used  in  the  study.  The 
broad,  dark  shaded  bands  on  the  images  are  a  result  of  the  technique  used 
to  photograph  the  screen  of  the  video  monitor  in  the  International  Imaging 
Systems  Model  70  image  processing  system  used  in  this  study.  Various 

textural  edges  in  these  images  are  used  in  the  studies  detailed  in  the 
following  sections. 


Figure  2.1  Seasat  Image  1 


Figure  2.2  Seasat  Image  2 


2.2  Erosion  Study 

Erosion  is  a  filtering  approach  applied  to  binary  images.  The  process 
is  summarized  in  Reference  6  as  follows.  The  basic  idea  is  to  define  a 
structural  element  as  a  set  of  resolution  cells  constituting  a  specific 
shape  such  as  a  line  or  square  and  to  generate  a  new  binary  image  by 
translating  the  structural  element  through  the  image  and  retaining  only 
those  pixels  as  l's  where  there  is  a  match  between  the  structural  element 
l's  and  the  image  l's.  The  process,  in  effect,  erodes  the  binary  image  as 
successive  translations  of  the  structural  element  are  made  through  the 
image.  This  process  is  shown  very  simply  in  Figure  2.3. 
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Eroded  by 


Figure  2.3  Erosion  Process 


The  textural  feature  obtained  for  each  translation  of  the  structural 
element  through  the  image  is  the  number  of  l's  left  after  each  cycle.  For 
binary  Images  this  is  the  same  as  the  area.  The  area  versus  the  number  of 
erosion  cycles  is  plotted  and  yields  what  is  called  the  covariance 
function. 

The  power  of  this  approach  is  that  it  emphasizes  the  shape  aspects  of 
the  tonal  primitives  of  the  texture.  It  has  found  wide  application  in  the 
analysis  of  microstructures.  Since  'the  texture  in  Figures  2.1  and  2.2  are 
very  fine,  i.e.,  the  tonal  primitives  are  very  small,  it  seemed  that  this 
approach  could  be  used  to  determine  shape  characteristics,  i.e., 
orientation,  width,  and  density  of  the  different  textural  edges.  The 
edges  considered  in  this  study  are  shown  in  Figure  2.4  and  Figure  2.5. 
The  coordinates  of  the  crosshairs  define  the  names  of  the  two  edges,  i.e., 
edge  (430,  312)  and  edge  (127,  122).  The  edges  were  partitioned  into  a 
32x32  image  for  the  analysis. 


Figure  2.4  Figure  2.5 
Location  of  Edge  (430,312)  Location  of  Edge  (127,122) 

The  important  parameters  in  the  erosion  procedure  are  the  binary  image 
and  the  shape  and  size  of  the  structural  element,  or  mask,  used  to  erode 
the  image.  The  binary  image  is  important  from  the  standpoint  of  the 
threshold  level  used  to  generate  it  from  the  original  grey  level  image. 
Figures  2.6,  2.7,  2.8  are  the  results  of  thresholding  Figure  2.1  at  levels 
where  10%,  20%,  and  30%,  respectively,  of  the  grey  levels  in  the  original 
are  above  the  threshold  level.  All  pixels  above  the  threshold  value  are 
set  to  1  in  the  binary  image  while  all  pixels  below  the  threshold  are  set 
to  0. 

The  structural  elements  used  in  this  study  were  primarily  chosen  to 
determine  how  well  the  orientation  and  density  of  the  edges  could  be 
determined.  Figures  2.9,  2.10,  2.11  show  the  effect  of  one,  two,  and 
three  erosion  cycles,  respectively,  on  Figure  2.6.  The  structural  element 
used  was  a  horizontal  line  three  elements  long,  i.e.,  [1,1,1].  The 
covariance  plots  resulting  from  the  complete  erosion  process  on  edge  (127, 
122)  and  edge  (430,  312)  in  Figures  2.6,  2.7,  and  2.8  are  shown  in  Figures 
2.12  and  2.13. 


8 


Figure  2.6 

Seasat  Inage  1  -  10%  Threshold 


Figure  2.7 

Seasat  Image  1  -  20%  Threshold 


Figure  2.8  Figure  2.9 

Seasat  Image  1  -  30%  Threshold  One  Erosion  -  10% .Threshold 


Figure  2.10 

Two  Erosions  -  10%  Threshold 


Figure  2.11 

Three  Erosions  -  10%  Threshold 


Figure  2.13  Covariance  for  Edge  (127,122) 


E^ge  (127,  122)  and  edge  (4.30,  312)  were  chosen  because  of  their 
decidedly  different  characteristics,  as  evident  in  Figures  2.4  and  2.5. 
An  examination  of  the  covariance  plots  in  Figures  2.12  and  2.13  shows  that 
erosion  with  the  three-wide,  horizontal  structural  element  results  in 
plots  that  are  markedly  different  a’so.  Edge  (430,  312)  is  oriented  at 
approximately  45  degrees  and  is  relatively  broad  hut  not  very  dense. 
Visually  it  is  not  perceived  as  having  any  horizontal  qualities  due, 
primarily,  to  its  density  and  orientation.  This  is  also  evident  on  the 
covariance  plot  of  Figure  2.13.  The  very  rapid  erosion  to  zero  in  five 
cycles  for  every  threshold  level  indicates  very  little  in  the  way  of 
horizontal  qualities.  The  large  area,  or  number  of  l's,  at  the  30% 
threshold  level  is  indicative  of  the  size  of  the  edge  hut,  due  to  its 
orientation,  the  horizontal  mask  quickly  shows  that  the  structural 
primitives  comprising  the  texture  of  this  edge  are  not  horizontal  in 
nature . 

Edge  (127,  122)  is  visually  nerceived  as  a  dense,  relatively  thin 
horizontal  edge.  Its  covariance  plot  in  Figure  2.12  immediately  reflects 
the  horizontal  orientation  through  the  large  number  of  cycles  necessary  to 
erode.  Even  more  important  is  the  piecewise  linear  nature  of  the  plots. 
Constant  slope  is  an  indication  that  the  same  number  of  pixels  are  being 
eroded  for  each  cycle.  This  implies  that  a  very  uniform  structure 
matching  the  orientation  of  the  structural  element  Is  present.  Figures 
2.*),  2.10,  and  2.11  also  show  that  during  the  erosion  process  the 
horizontal  qualities  of  the  edge  were  amplified  by  the  horizontal  mask. 

Two  other  structural  elements  masks  were  applied  to  edge  (430,  312). 
These  elements  were  oriented  at  135  and  45  degrees.  The  covariance  plots 
for  them  are  shown  in  Figure  2,14  and  Figure  2.15.  The  same  properties 
are  evident  as  discussed  previously  for  the  horizontal  mask.  The  45 
degree  mask  Immediately  shows  the  orientation  of  the  edge.  The  piecewise 
linear  structure  emphasizes  the  orientation  match.  During  the  erosion 
process  the  45  degree  property  was  emphasized  by  the  mask. 
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Figure  2.14  Erosion  Cvcles  Covariance  for  Edge  (430,  312),  -  Mask 
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Figure  2.15  Erosion  Cvcles  Covariance  for  Edge  (430,  312),  -  Mask 
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2.2.1  Conclusions 


The  eros'on  process  was  very  good  for  emphasising  the  hasic  properties 
of  the  edges  considered  in  this  study.  Emphasis  was  placed  on  orientation 
and  density  properties  in  this  study.  There  Is  no  reason  that  periodicity 
of  texture  structures  cannot  he  Identified  by  using  structural  elements 
designed  to  accent  that  property.  Such  masks  would  consist  of  l's 
separated  by  0's  with  the  blank  space  determining  the  period. 

The  texture  images  resulting  from  this  SAR  sensor  seem  to  lend 
themselves  very  well  to  the  texture  analysis  approach.  One  problem  area 
is  that  it  Is  an  Iterative  process  requiring  sometimes  many  passes.  One 
possible  very  good  use  for  this  process  would  be  to  use  It  as  a 
preprocessing  step.  A  few  erosion  cvcles  could  be  performed  to  accent  any 
areas  of  an  image  that  match  the  properties  the  structural  element  was 
designed  to  reveal.  Then  the  resulting  Image  could  be  used  for  further 
processes  but  the  properties  of  interest  would  now  dominate  other 
qualities  or  features  in  the  image.  This  could  aid  the  feature  selection 
process  for  representing  image  data. 

2.3  Statistics  Study 

8 

K.I.  Laws  developed  and  reported  on  "texture  energy"  transforms 
which  performed  better  than  co-occurance  statistical  approaches.  For  a 
zero  mean  field,  the  texture  energy  measure  is  the  standard  deviation 
since  the  variance  would  be  the  average  of  squared  signal  values,  an 
energy  measure  in  the  formal  sense  of  the  word.  If  the  Image  had  been 
previously  filtered,  the  texture  energy  measures  the  local  energy  within 
the  pass  band. 

Since  either  the  variance  or  the  standard  deviation  alone  has  been 
shown  to  be  sufficient  to  extract  texture  information,  a  statistics1  study 
of  two  types  of  textural  edges  was  performed.  The  study  consisted  of 
determining  the  mean  and  standard  deviation  of  a  local  area  as  that  area 
was  moved  across  the  textural  edges.  The  edges  used  in  this  study  are 
shown  In  Figure  2.2  and  Figure  2.16. 
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Since  the  mean  and  standard  deviation  are  local  operations,  the  size 
of  the  local  area,  or  window,  is  an  important  parameter.  Two  sizes  of 
windows  were  used,  a  32x32  and  a  16x16  area.  There  was  no  particular 
reason  for  choosing  these  sizes  other  than  they  fit  within  the  texturallv 
distinct  areas  that  comprised  the  region  around  the  edges.  For  faster 
iterative  operation,  smaller  sizes  would  be  more  appropriate. 

The  location  of  the  window  is  defined  by  the  upper  left  hand  corner 
pixel  co-ordinates.  When  the  window  was  moved  through  the  edges  it  was 
placed  to  the  left  of  the  edges  and  moved  completely  through  them.  This 
is  normally  done  in  a  raster-type  scan  in  most  hardware.  For  an 
Interesting  excursion,  a  path  perpendicular  to  the  orientation  of  the  edge 
was  used  to  determine  if  being  able  to  depart  from  the  normal  raster-scan 
method  of  convolving  a  local  operator  through  an  image  could  be  beneficial. 

Statistics  were  also  determined  for  the  edges  after  they  had  been 
filtered  with  a  Sobel  operator.  The  Sobel  gradient  is  an  edge  detection 
operation.  It  has  been  used  in  texture  discrimination  studies  with  good 
results.  It  1 s  a  nonlinear  3x3  operator  which  is  defined  by  the  following 
masks: 
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For  each  pixel  the  Sobel  magnitude  is  determined  by 


In  this  study,  the  square  of  the  Sobel  magnitude  was  used  and  the 
operation  was  labeled  SOBELSQ.  Figures  2.17  and  2.18  show  the  results  of 
applying  this  mask  to  the  edges  of  interest. 

Figures  2.18  -  2.24  show  the  results  of  the  statistical  measurements. 
As  is  evident  in  Figure  2.19,  the  32x32  mask  size  does  not  descrtminate 
the  narrow  edge  (430,312).  This  is  a  result  of  the  edge  comprising  too 
small  a  percentage  of  the  mask  area  as  the  mask  is  moved  through  it. 
Figure  2.20  shows  that  the  smaller  mask  size  greatly  improves  the  edge 
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descriminat ion  capability.  It  also  shows  that  the  SOBELSQ  operation 
provides  little  improvement.  Figure  2.21  indicates  that  an  appreciable 
improvement  in  locating  the  edge  boundary  can  be  obtained  by  moving  the 
mask  on  a  path  perpendicular  to  the  edge.  This  results  from  the  increased 
percentage  the  edge  has  in  the  mask  as  it  first  enters  the  mask. 

Figures  2.22  -  2.24  show  the  results  on  edge  (72,372).  This  edge  is 
different  from  the  narrow  edge  (430,312).  It  has  more  width  and, 
therefore,  should  show  a  double  mode  characteristic  as  the  mask  is  passed 
through  it.  In  fact,  both  the  32x32  and  16x16  masks  do  not  meet  this 
expectation,  as  shown  in  Figures  2.22  and  2.23.  The  standard  deviation 
measure  is  enhanced  by  the  Sobel  operation  in  the  16x16  case.  Figure  2.24 
shows  the  double  mode  characteristic  is  obtained  when  the  mask  is  moved  on 
a  perpendicular  path  to  the  edge  orientation. 

2.3.1  Conclusions 

As  expected,  mask  size  and  path  direction  are  important  in  determining 
the  effectiveness  of  a  filtering  operation.  The  usefulness  of  the 
variance  and  standard  deviation  as  texture  measures  for  these  images  is 
not  completely  evident  from  the  limited  results.  Certainly  under  the 
right  conditions  of  mask  size  and  filter  preprocessing,  such  as  a  Sobel 
magnitude  operation,  the  usefulness  of  the  variance  as  a  texture  measure 
could  be  enhanced.  Certainly  the  present  results  do  not  conclusively 
point  to  the  variance  aa  the  sought  for  texture  information  indicator  for 
these  texture  images.  Without  a  priori  knowledge  of  the  texture 
characteristics,  however,  it  remains  an  effective  texture  measure. 
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Figure  2.21  16x16  Statistics  for  Perpendicular  Path 
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Figure  2.22  32x32  Statistics  for  Edge  (72,  372) 
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SECTION  III.  REACHABLE-SET  GUIDANCE  TECHNIQUE 


3.0  Introduction 

Guidance  of  short  range  air-launched  missiles  has  been  based  largely 

on  a  static  technology  for  the  last  decade.  Conventional  guidance  of 

"fire  and  forget"  weapons  has  been  dominated  by  variations  of  proportional 

and  pursuit  navigation,  founded  in  optimal  control  for  non-maneuvering 
9-13 

Interception  .  They  require  simple  angle  measurements,  and  are 
easily  Implemented  In  analog  hardware.  Within  terminal  saturation 
constraints  such  techniques  provide  adequate  accuracy  for  scenarios 
involving  non-maneuvering  vehicles.  Yet,  an  intelligent  target  capable  of 
radical  maneuvers  can  deceive  a  conventional  guidance  strategy  bv 
purposely  Inducing  terminal  saturation  in  the  final  seconds  of  the 
encounter. 

In  the  presence  of  maneuvers,  formulation  of  an  effective  guidance 

strategy  becomes  far  more  complex^  Control  effort  must  in  some 

sense  anticipate  target  trajectories  by  modeling  maneuver  capabilities  and 

the  target  response,  both  deterministic  and  random,  to  the  closing 

missile.  Stated  in  this  framework,  the  problem  reduces  to  a  differential 
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games  formulation,  defining  optimal  evasion/pursuit  strategies  .  Yet, 

the  solution,  even  for  very  contrived  scenarios,  leads  to  a  significant 
numerical  burden,  unsuitable  for  practical  usage. 

While  such  an  elaborate  approach  will  be  of  doubtful  value  in  this 

context  for  some  time  to  come,  it  does  serve  to  Indicate  that  onboard 
Intelligence  can  be  used  to  greatly  improve  a  missile's  advantage  over  Its 
adversary.  Indeed,  the  capabilities  of  digital  hardware  have  matured  to 
the  point  where  serious  consideration  must  be  given  to  advances  in 
guidance  strategy  over  the  conventional  techniques.  Of  particular  value 
would  be  a  study  of  the  tradeoffs  possible  between  level  of  intelligence 
(l.e.,  hardware  capabilities)  and  overall  missile  performance  (i.e.,  miss 
distance  and  aspect  angle).  An  extensive  evaluation  of  this  nature  would 
illuminate  key  factors  necessary  for  the  development  of  future  weapons 
delivery  systems.  By  reducing  computational  requirements  to  a  common 


20 


denominator,  t.e.,  currently  available  hardware,  practicality  with  respect 

to  physical  size  limitations  can  be  assessed;  this  also  allows 

extrapolation  of  practicality  into  the  future  with  projected  advances  in 

microelectronics.  Of  course,  it  is  the  performance  of  any  given  control 

strategy  that  ultimately  determines  if  the  necessary  hardware  is  warranted. 

The  work  effort  described  by  this  report  does  not  attempt  such  a 

survey,  but  rather,  a  small  fraction  of  such  a  study  was  conducted;  a 

single  promising  guidance  technique  was  examined  and  compared  with 

benchmark  simulation  tests  of  conventional  guidance  in  several  encounter 

scenarios.  Then,  using  currently  available  technology,  the  architecture 

of  the  required  hardware  was  examined.  The  results  indicate  that  using 

current  technology,  substantial  improvements  in  missile  performance  can  be 

realized.  Of  course,  this  is  simply  a  single  point  of  the  overall  study, 

and  does  not  pretend  to  proclaim  the  best  currently  available  technique. 

The  subject  of  this  evaluation  is  a  guidance  law  developed  in 

Reference  20,  based  on  the  concept  of  reachable  set  theory  for  dynamic 
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systems  .  It  was  chosen  because  it  is  representative  of  a  class  of 

guidance  strategies  that  could  be  of  immediate  value,  i.e.,  it  (1)  has  an 

intuitive  structure,  (2)  calls  for  moderate  computation  in  the  form  of  a 

systematic  search,  and  (3)  is  suited  for  a  sampled  data  context. 

Variations  can  be  appended  to  the  basic  law  to  adapt  it  to  other  types  of 
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encounters  by  using  a  cost  function  based  on  physical  limitations 

The  following  sections  present  details  regarding:  Encounter  Model 
(Section  3.1),  Conventional  Guidance  Implementation  (Section  3.2), 
Advanced  Guidance  Implementation  (Sections  3.3  through  3.5),  and 
Observations  of  Performance  From  Simulations  (Section  3.6). 

3.1  Encounter  Model 

Before  discussing  the  guidance  techniques  in  detail,  it  is  necessary 
to  describe  the  model  of  the  target  and  missile  behavior.  For  the  most 
part,  the  assumptions  and  numerical  values  used  in  reference  20  were  drawn 
upon.  The  target  was  allowed  a  constant  forward  speed  of  1000  ft/sec; 
maneuver-induced  drag  was  assumed  compensated  by  forward  thrust.  The 
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maximum  turn  was  governed  by  a  6  G  constraint  on  normal  acceleration.  Due 
to  the  short  duration  of  the  close-range  encounter,  the  target  was 
constrained  to  maneuver  in  a  single  plane  having  a  "tilt"  angle  6  with 
respect  to  its  initial  velocity  vector,  as  shown  in  Figure  3.1. 


Figure  3.1  Maneuver  Plane  Definition 


The  missile  model  was  somewhat 
an  initial  speed  of  1000  ft/sec, 
thrust  was  given  as  4700  lb,  with 
was  50  lb  of  the  initial  165  lb. 
vector,  magnitude  ax  ,  at  an  angle 
shown  in  Figure  3.2. 


more  complex.  It  was  assumed  to  have 
launched  from  its  host  aircraft.  The 
a  burn  tiem  of  2.6  sec.  The  fuel  load 
Guidance  was  by  means  of  a  normal  lift 
0  with  respect  to  the  missile  body  as 
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Figure  3.2  Missile  Acceleration  Vector 
In  the  interest  of  practicality,  a  parabolic  drag  force  law  was  used 
0  -  h  Ado  +  K2  CL 

V 


where 

S)0  *  2.3  -  zero  lift  drag  coefficient 
*  .0025  =*  induced  drag  coefficient 
”  proportionality  factor  =  .001 

■  proportionality  factor  =  1000 

w  ■  missile  weight,  function  of  time 

g  ■  gravitational  acceleration 

■  missile  speed 
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Scenarios  were  defined  by  the  relative  positions  of  the  target  and 
missile  and  their  respective  velocity  vectors.  The  target  was  then 
allowed  a  maneuver  strategy  bounded  by  6  G's  in  turn  rate,  at  a  fixed 
angle  of  tilt.  Missile  guidance  was  done  by  choice  of  lift  vector  as  a 
function  of  closure. 

3 . 2  Implementation  of  Conventional  Guidance 

For  the  purposes  of  simulation  and  benchmark  evaluation,  conventional 
guidance  was  represented  bv  two  techniques;  proportional  navigation  and 
pursuit  guidance.  Linear  combinations  of  their  respective  components  can 
be  used  where  the  weighting  constants  may  be  functions  of  range,  i.e., 
favoring  pursuit  initially  and  proportional  on  final  approach^  .  Here, 
each  was  used  in  its  purest  form. 


3.2.1  Proportional 

The  magnitude  of  the  normal  acceleration  for  proportional  guidance  is 
given  as 

<c)  -  c  IW^J 

where 

# 

®L0S^  =  rotational  rate,  with  respect  to  inertial  space,  of 
the  line  of  sight  angle  to  the  target 
Vc  =  closing  speed 

C  =  navigation  constant  (normally  between  3  and  6) 


Missile 


I  Inertial 


Figure  3.3  Proportional  Guidance 


In  three  dimensions,  the  angle  of  the  missile  acceleration  orientation 
i 8  chosen  to  rotate  the  missile  velocity  vector  toward  the  target's 
relative  displacement  vector. 

While  in  practice  this  calculation  is  done  in  analog  hardware,  for 
simulation  purposes  the  necessary  derivative  was  approximated  by  a 
sampled-data  finite  difference. 
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where  T  was  a  small  sampling  interval.  At  each  such  iteration,  a  new 
control  effort  was  computed  and  used  to  drive  the  system's  dynamic 
equations. 


3.2.2  Pursuit 


For  this  second  case  the  magnitude  of  control  is  given  by 


*X<t>  ■  Cl0LOOK(t)lVc 


with  orientation  0  chosen  as  in  the  previous  case.  The  angle 
measured  from  the  missile  body  axis  to  the  target  location. 


LOOK 


is 


Figure  3.4  Pursuit  Guidance 

Again,  this  angle  was  sampled  at  uniform  intervals  for  the  purposes  of 
simulation. 

As  mentioned  earlier,  such  techniques  have  long  been  used  with 
success.  They  prove  adequate  for  launchings  where  terminal  saturation  is 
avoided  due  to  sluggish  target  evasion  or  close  range.  The  principal 
advantages  lie  in  its  means  of  implementation;  the  measurements  required 
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are  simple,  i.e.,  only  a  reasonable  guess  at  closing  speed  and  an  accurate 
estimate  of  the  displacement  angle.  The  simple  measurements  and  their 
associated  sensors,  coupled  with  the  analog  control  hardware,  make  such 
navigation  schemes  attractive  from  the  point  of  view  of  cost  and  physical 
size. 

Yet,  because  these  schemes  are  based  on  non-evasive  targets,  they  tend 
to  de-emphastze  any  initial  launch  advantage,  and  prefer  to  postpone 
offensive  counter-maneuvering  until  late  in  the  scenario.  That  is,  since 
the  controller  does  not  expect  changes  in  the  target  trajectory,  it 
responds  only  after  a  maneuver  becomes  evident  at  the  angle  sensor.  This 
is  often  too  late  for  adequate  course  correction,  and  invariably  results 
in  terminal  saturation  and  miss  distance  dependent  on  the  agility  of  the 
target.  One  means  of  dealing  with  this  terminal  miss  effect  is  to 
increase  the  warhead  size  and  the  kill  radius. 

3.3  Advanced  Guidance 


The  terminal  effects  associated  with  conventional  guidance  are  largely 
responsible  for  the  interest  in  more  sophisticated  navigation  techniques 
capable  of  anticipating  and  responding  to  target  maneuvers.  A  missile 
launched  with  a  high  kill  probability  means  that  interception  is  highly 
likely  for  all  valid  target  maneuvers.  That  is,  the  target's  position  in 
space  and  time  will  be  contained  in  the  set  of  all  points  reachable  by  the 
missile;  as  time-to-go  decreases,  this  reachable  set  shrinks.  The 
guidance  controller  must  maintain  the  initial  advantage  by  anticipating 
target  attempts  to  exit  the  missile's  shrinking  reachable  set.  For  this, 
the  controller  should  examine  all  valid  target  trajectories  and  respond  as 
if  the  worst  one  (from  the  missile's  viewpoint)  were  to  be  used. 

Clearly,  this  is  an  ill-posed  problem  which  can  be  rendered  tractable 
by  the  quantitative  observation  in  reference  20.  Using  a  differential 
games  analysis  of  this  framework  of  assumptions,  it  can  be  shown  that  the 
target  maximizes  time-to-go  when  its  maneuver  is  a  maximal  acceleration 
turn.  For  our  purposes,  the  tilt  angle  6  of  this  maneuver  is  unimportant; 
a  complicated  function  of  relative  position  and  attitude.  Using  this 
observation,  a  tractable  guidance  scheme  can  then  be  implemented.  In 
brief,  the  controller  determines  the  target's  trajectory,  assuming  the 
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target  makes  Its  maximum  turn  at  some  selected  tilt  angle  5  •  Upon 
examination  of  a  set  of  such  angles,  the  missile  responds  as  if  the  target 
were  to  choose  the  "worst  case",  i.e.,  anticipating  the  optimal  maneuver 
for  the  target.  After  a  brief  interval,  e.g.,  determined  by  computational 
requirements,  the  process  is  repeated  using  updated  position  and  attitude 
information,  yielding  a  "closed-loop”  implementation. 

The  detailed  operation,  including  mathematical  particulars,  is  as 
follows : 

Given  observations  of  relative  displacement,  missile  and  target 
heading  and  speed,  assuming  a  relative  coordinate  system  to  be  described 
later,  then 

a.  the  controller  assumes  that  the  target  chooses  a  maneuver 
plane  defined  by  5  and  uses  its  optimal  maneuvers  to  maximize  time-to-go. 

b.  Next,  the  necessary  missile  heading  coordinates  (  <f>  ,  0  )  to 
effect  an  intercept  are  computed.  This  is  done  using  the  following  time 
function  expressions  for  relative  Cartesian  displacement: 

Ax(t)  =  Ax(0)  +  V  t  cos(4>)  cos(0)  -  -  sin(V  t  /a) 
m  at  t  t 

Ay(t)  =  Ay(0)  +  V  t  sin(<f))  cos (0)  +  —  cos(6)  [cos(v  t  /a  )  -  1] 
m  at  t  t 

Az(t)  =  Az(0)  +  V  t  sin(<())  +  —  sin(6)  [cos(v  t  /a  )  -  1] 
m  a_  t  t 


where 

V  *  target  speed 

a  ”  target  turn  rate  (maximum! 

V  *  missile  speed 
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The  closed  form  Integration  results  by  assuming  constants  for  speed, 
maneuver  angle,  and  heading  angle.  At  the  interception  time  t^ 

Ax  -  Ay  •  Az  ”  0 

By  defining  A  ■  V  tr  as  the  distance  traveled  by  the  missile,  both  <p 

m  r 

and  8  can  be  eliminated  to  yield  a  single  scalar  equation  in  the 
tlme-to-go: 

f(tj  »  Jl2  -  [Ax(0)~  —  sin  (V.t  /a.)]2 

£  at  t  i  t 

-  [Ay(0)~  j  cos(6)  (cos(Vttf /at)-l) ]2 

-  [Az(0)~  —  sin(6)  (cos(V  t  /a  )-l)]2 

at  t  r  t 

This  can  be  solved  for  the  positive  root  t^  using  a  Newton-Raphson 
search;  then  the  actual  missile  heading  can  be  found  using  closed-form 
evaluations : 

♦  “  W  =  8in 


[A  z(0)+  f 


sin  (6) 


[cos(Vttf/at) 


-  D] 


v 

m  f 


0  ”  f8^tf^  =  COS 


At  this  point,  the  velocity  vector  that  the  missile  should  have  for 

interception  is  known,  (V  <j>  9  ),  given  the  specific  target  maneuver 

tn 

angle  6 . 


[Ax(0)+  —  sin  (V  t  /a  )] 
at  t  r  r 

V  t.  cos  ( 4>) 
m  f 


c.  Given  the  actual  missile  heading,  the  acceleration  necessary 

to  yield  an  average  velocity  vector  of  (V  <p  ,0)  over  the  interval 
- * —  m 

(0,  tj.)  is  approximated  by 


a j_  (6) 


where  §  is  the  angle  separating  initial  and  desired  velocity  vectors. 
This  relationship  is  depicted  in  Figure  3.5  and  represents  the  control 
effort  necessary  to  respond  to  a  specific  maneuver. 


Necessary 
Change  in 
Velocity 
Av 


Figure  3.5  Required  Missile  Velocity  Change 

d.  Conceptually,  a^functlon  a  (<$ )  exists,  giving  the  necessary 
missile  control  effort  as  a  function  of  target  maneuver  angle  6  .  The 
maximum  of  this  function  defines  the  worst  case  evasion  that  the  target 
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can  choose.  Hence,  the  missile  controller  anticipates  this  as  the 

★ 

maneuver,  and  responds  with  normal  acceleration  a^ ,  oriented  at  the 
angle  of  vector  AV  defined  above.  The  computation  involved  will  be 
detailed  in  the  next  section. 


-TT  6*  IT 

Maneuver  Angle  5 


Figure  3.6  Missile  Acceleration  Function 

Note  that  an  interesting  feature  emerges  for  such  an  analysis.  At  any 
given  instant,  the  function  a±  (6)  summarizes  the  advantage  that  the 
missile  has  over  its  target.  That  is,  if  a  means  of  successful  evasion 
exists,  then  for  some  angle  6,  aj  (6  )  exceeds  the  maximum  allowable  normal 
acceleration  of  the  missile.  Such  information  could  be  particularly 
valuable  to  a  pilot  if  presented  as  a  "probability  of  hit"  measure. 

While  computation  and  architecture  will  be  discussed  in  subsequent 
sections,  one  consideration  must  be  mentioned  here.  This  guidance  scheme 
in  its  closed-loop  form  must  actually  be  updated  in  a  continuous  fashion. 
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Yet  due  to  computational  operations,  it  becomes  sampled  data  control,  with 
the  update  or  sampling  interval  determined  by  the  necessary  computation. 


In  this  section,  the  requirements  of  the  advanced  reachable-set  based 
guidance  scheme  are  discussed.  First,  an  overall  description  of  the 
software  is  presented,  addressing  the  structure  of  the  Fortran  listing  of 
GUIDE  found  in  Appendix  A.  Then,  using  findings  from  tests  using  this 
software,  currently  available  hardware  is  evaluated. 

A  simplified  flowchart  in  Figure  3.7  is  included  here  to  aid  in  the 
description  of  the  Fortran  listing  in  the  Appendix.  Upon  first  entry, 
initialization  steps  are  encountered,  allowing  the  user  to  set 
interactively  a  number  of  options  and  parameters.  At  run  time,  this 
section  (lines  34-74)  is  skipped,  and  actual  control  computation  proceeds 
as  follows: 


1)  First,  the  elapsed  time  from  the  previous  control  calculation 
is  checked  against  the  specified  update  interval.  If  insufficient  time 
has  passed,  there  is  an  Immediate  return  hack  to  the  calling  program.  The 
control  effort  previously  determined  is  maintained. 

2)  When  the  update  interval  has  elapsed,  a  new  control  effort  is 
determined  from  current  measurements.  Passed  to  the  routine  are  the 
actual  states  of  missile  and  target  with  respect  to  an  inertial 
reference.  The  pairs  of  six-dimensional  vectors  specifically  contain  (X, 
Y,  Z)  position  values,  followed  by  spherical  velocity  information  (V,<J)  ,0), 
i.e.,  speed,  azimuth  and  elevation  of  the  velocity  vector.  These  are 
transformed  by  linear  algebraic  rotation  to  a  missile-centered  coordinate 
system.  Use  of  this  coordinate  system  as  a  basis  for  necessary 
measurements  eliminates  the  need  for  a  strapped-down  Inertial  reference  in 
the  missile.  Assuming  the  turn  rates  are  slow  with  respect  to  the  update 
interval,  only  moderate  degradation  in  the  overall  performance  is 
experienced . 

The  new  coordinate  system  is  defined  by  an  X-axis  in  the  direction  of 
the  missile  axis,  and  its  positive  Y-axis  lying  in  the  plane  of  the  target 
point . 


3)  The  five  measurements  in  this  system  (range,  target  azimuth, 
target  speed,  velocity  vector  azimuth  and  elevation)  are  appropriately 
corrupted  by  random  measurement  noise. 

At  this  point  the  target  measurements  have  been  conditioned  as  if  the 
missile  sensors  had  gathered  them;  i.e.,  they  represent  the  information 
available  from  sensors  having  only  the  missile  and  target  as  directional 
reference.  Thus,  the  code  in  the  subroutine  to  this  point  (line  100)  is 
simply  overhead  computation.  Actual  control  computation  made  bv  the 
missile  begins  at  line  130.  Certain  sections  of  the  code  associated  with 


the  missile 

are 

also 

overhead,  performing  certain 

initialization 

computations. 

As 

such, 

they 

are  executed 

only  one  time 

per  pass  so 

efficient  coding 

was  not 

felt 

necessary. 

The  computation 

bound  loops, 

however,  must 

he 

studied 

for 

improvement 

in  efficiency 

before  actual 

implementation. 
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4)  To  begin,  the  missile-centered  measurements  were  transformed 
to  a  target-centered  system,  defined  analogously  to  the  missile-centered 
system.  This  step  is  necessary  due  to  the  problem  formulation  described 
in  Section  3.3.  This  is  done  in  lines  130  through  165. 


5)  At  this  point,  the  search  for  the  worst  case  control  effort 

★ 

begins.  A  window  is  defined  straddling  the  angle  6  determined  as  the 
target's  worst  maneuver  at  the  previous  update  time.  In  practice, 
assuming  a  sufficiently  fast  update  rate,  the  tilt  angle  of  worst  maneuver 
changes  slowly,  so  will  remain  within  such  an  interval.  (For  the 
simulations  conducted,  an  interval  of  30°  was  used.)  For  the  two 
extreme  maneuver  angles,  <$i  =  6  -  -2»^2=6  +  y  the  solution  is  found  for 

the  necessary  missile  velocity  vector  direction  to  effect  an  intercept. 
As  described  earlier,  this  is  done  using  a  Newton-Raphson  search  for  the 
time-to-go,  t-,  and  then  solving  a  closed-form  expression  for  <)>  and  0  . 
Assuming  good  starting  values  (the  previous  value  found  for  t^  is 
adequate)  convergence  will  require  only  two  or  three  Iterations.  For  each 
window  edge,  the  necessary  control  effort  a^  is  found,  again  described  in 
Section  3.3. 


6)  The  process  of  search  and  evaluation  is  repeated  using  the 
window's  midpoint. 

7)  The  smaller  of  the  control  efforts  computed  for  the  window 
edges  is  determined  and  the  corresponding  tilt  angle  abandoned  in  favor  of 
the  midpoint  angle.  That  is,  the  window  is  halved  by  rejecting  the 
maneuver  angle  that  represents  the  lesser  threat.  The  control  effort  as 
well  as  the  desired  velocity  vector  parameters  are  stored  for  the  new 
window  edge. 


8)  The  new  window  width  is  checked  against  a  prespecified 
threshold.  If  it  exceeds  the  threshold,  the  binary  search  continues  by 
repeating  step  (6)  above.  If  it  is  Indqed  less,  the  maximization  of 
control  effort  is  complete. 


9)  At  this  point,  all  pertinent  Information  about  the  worst 
possible  maneuver  is  available.  Specifically,  this  includes  the  time  to 
intercept,  t^,  the  required  control  effort,  aj_  ,  and  the  desired  missile 
velocity  vector,  (V,  (p  ,  0  ) .  The  vector  is  transformed  back  to 
missile-centered  system. 

10)  Finally,  the  angle  of  application  for  the  control  vector 

is  determined.  At  this  point  the  missile  has  sufficient  information  to 

determine  the  necessary  control  surface  deflections  to  generate  its 
acceleration  vector. 

11)  The  last  section  of  the  subroutine  (lines  297  through  309) 
simply  map  the  acceleration  vector  back  to  the  inertial  system,  for 
compatibility  with  the  calling  program. 

3.5  Computational  Requirements 

The  proposed  algorithm  is  heavily  arithmetic  bound.  As  a  result, 
computational  requirements  are  easily  estimated  since  the  time  required 
for  floating-point  arithmetic  will  be  predominant  and  will  be  a  good 
estimate  of  total  computation  time  required.  Additionally,  software 
emulated  floating-point  arithmetic  would  obviously  not  be  acceptable. 

The  algorithm  was  divided  into  its  major  parts  as  shown  in  the  block 
diagram  of  Figure  3.8  and  computational  times  estimated  for  each  part.  An 
8086  microcomputer  operating  at  5  MHz  with  an  auxiliary  8087 
floating-point  processor  was  used  for  all  time  estimates.  Estimates  were 
made  by  counting  floating-point  operations  (including  load  and  store)  in 
the  original  Fortran  program  and  multiplying  by  the  appropriate  8086/87 
instruction  time.  Pessimistic  estimates  were  made  at  all  times. 

In  Figure  3.8,  the  start-up  search  was  ignored  since  it  is  only  made 
once  and  contributes  little  to  the  com""' ational  problem's  dynamics. 
Times  for  each  of  the  subprocesses  are  shown  in  Figure  3.8.  The  total 
computational  time  required  by  the  8086/87  processor  is  given  by: 

37,877  +  14,586  N  +  9,304  M  +  7,293  M  N 
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where  N  is  the  number  of  Newton  iterations  needed  to  calculate  the 
direction  of  the  velocity  vector  for  intercept  and  M  is  the  number  of 
binary  search  passes  required  to  search  for  the  maneuver  angle.  Typical 
values  are  three  Newton  iterations  and  M“4  corresponding  to  dividing  the 
maneuver  angle  search  window  into  16  segments.  For  these  values,  the 

total  computational  time  is  approximately  206.4  milliseconds,  about  an 
order  of  magnitude  slower  than  desired. 

3.5.1  Trigonometric  Look-Up  Table 

The  algorithm  was  searched  for  significant  opportunities  for 

improvement.  Although  the  8087  floating-point  processor  was  used  for 
calculating  a  variety  of  functions,  the  sine  and  cosine  functions  were 
especially  prevalent  and  time  consuming.  The  8087  calculates  these  two 
functions  from  the  tangent  function  via  trigonometric  identities.  Another 
time  estimate  was  calculated  using  the  8087  for  all  calculations  except 
for  9ine  and  cosine  evaluations.  These  latter  functions  were  assumed  to 
be  stored  in  a  look-up  table  in  read-only-memory.  The  results  of  this 
estimate  are  also  shown  in  Figure  3.8.  The  total  computation  time 

required  using  the  8086/87  and  a  look-up  table  for  sine  and  cosine  is 

given  by: 


17,893  +  5,766  N  +  4,024  M  +  2,883  M  N 

For  the  same  conditions,  N=3  and  M=4,  as  in  the  previous  example,  total 
computation  time  was  approximately  85.9  milliseconds.  Although  a 
significant  improvement  over  the  original  estimate  of  206.4  milliseconds 
without  the  look-up  table,  further  reduction  of  computational  time  was 
desirable  to  improve  the  performance  of  the  algorithm. 

3.5.2  Multiple  Processors 

The  next  logical  step  was  to  attempt  configurations  of  multiple 
8086/87  processors.  Simulation  showed  that  end-to-end  computational  delay 
was  critical  and  not  sampling  rate.  Thus,  pipelined  processor 
configurations  were  ruled  out  as  not  reducing  end-to-end  delay  but  only 


increasing  sampling  rate.  Parallel  computation  was  clearly  necessary. 
Examining  Figure  3.8,  most  computation  Is  clearly  required  In  setting  up 
the  search  window  and  searching  for  the  maneuver  angle.  Fortunately,  each 
of  these  two  tasks  could  he  configured  to  allow  parallel  processing.  Each 
edge  of  the  search  window  could  be  found  independently.  Each  half  of  the 
search  window  could  be  searched  independently.  Figure  3.9  is  the  block 
diagram  for  the  Implementation  of  the  algorithm.  The  time  estimates  in 
Figure  3.9  assume  two  completely  independent  8088/87  systems,  each  with 
its  own  sine  and  cosine  look-up  table.  These  systems  are  loosely 
coupled.  The  total  computational  time  for  the  dual  processor  system  is 
given  by: 

12,554  +  2,883  N  +  2,012  M  +  1,442  M  N 

Again,  for  N=3  and  M=»4,  the  total  computation  time  is  approximately  46.6 
milliseconds,  a  reasonable  performance  for  this  algorithm.  As  an  added 
bonus,  the  second  processor  could  be  used  for  sensor  and  actuator 
conditioning  as  shown  in  Figure  3.10.  The  times  when  the  second  processor 
would  not  be  needed  by  the  control  algorithm  are  exactly  these  times  when 
data  are  input  and  output. 

3.5.3  Future  Refinements 

Clearly,  all  three  performance  estimates  are  encouraging.  Even  the 
single  processor  estimate  of  206.4  milliseconds  is  sufficiently  fast  that 
a  newer-generatlon  processor  will  be  able  to  reduce  this  time  to  an 
acceptable  value.  Cost  and  performance  figures,  normalized  to  a  single 
processor  system  without  look-up  table,  are  shown  in  Table  3.1.  Notice 
that  although  the  dual  processor  system's  performance  is  greater,  its  cost 
performance  ratio  Is  actually  greater  than  that  of  the  single  processor 
with  look-up  table.  A  very  desirable  investigation  would  be  to  actually 
code  and  test  this  algorithm  on  an  8086  with  8087  floating-point  processor 
and  look-up  table  for  sine  and  cosine.  It  is  likely  that  an  improvement 
of  a  factor  of  2  over  the  pessimistic  estimates  could  be  found.  If  the 
performance  of  the  single  processor  system  could  be  so  improved,  the 
difficulties  and  costs  of  the  dual  processor  could  be  avoided.  In  anv 


case,  the  actual  Implementation  of  the  single  processor  system  would 
provide  more  accurate  data  on  which  to  base  performance  estimates  of  other 
systems. 


Figure  3.8  Single  Processor  Solutions 


Table  3.1.  Cost  and  Performance  of  Architectural  Alternatives 


SYSTEM 

PERFORMANCE 

COST 

COST/ PERFORMANCE 

8086/87  Processor 

1.0 

1.0 

1.0 

8086/87  Processor 
with  Look-up  Table 

2.4 

1.1 

0.46 

DUAL  8086/87  Processors 
with  Look-up  Table 

4.4 

2.2 

0.5 

3.6  Performance 

The  achievable  performance  of  an  advanced  guidance  technique 
determines  if  the  cost  of  additional  sensors  and  computational  hardware  Is 
warranted.  For  the  sake  of  comparison  with  conventional  guidance,  several 
classes  of  tests  were  simulated  using  the  target  and  missile  models 
described  earlier.  In  this  way,  the  robustness  of  the  reachable  set 
approach  could  be  assessed  by  systematically  degrading  the  assumptions  and 
measurements  entering  into  its  formulation.  Two  representative  encounters 
were  used,  shown  in  Figure  3.11:  Scenario  two  Involved  a  rear  attack  with 
the  target  maneuvering  by  pulling  up  at  two  o'clock,  l.e.,  a  low  crossing 
rate.  Scenario  five  dealt  with  a  side  attack  on  a  target  maneuvering  as 
before.  Involving  a  high  crossing  rate. 
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3.6.1  Perfect  Measurements 

To  evaluate  the  potential  for  success,  simulations  were  first  run 
using  perfect  measurement  Information.  The  resulting  miss  distances  are 
given  in  Table  3.2.  While  insignificant  computational  delay  was  assumed, 
a  sampling  or  update  rate  of  100  msec  was  used  for  the  reachable  set 
approach.  It  can  be  seen,  as  would  be  expected,  that  performance  of 
conventional  techniques  is  severely  degraded  by  high  crossing  rate.  In 
both  cases,  the  terminal  efforts  required  by  conventional  guidance 
saturated  the  allowable  limits  of  the  missile.  On  the  other  hand,  the 
advanced  approach  gave  a  control  effort  well  distributed  over  the 
encounter's  duration,  indicating  a  high  degree  of  anticipation,  as  noted 
in  reference  20.  Ideally,  the  effort  should  be  monotonically 

non-increasing,  whereas  in  practice  a  small  degree  of  "upturning"  of 
control  effort  is  witnessed  on  final  approach.  In  both  scenarios  the 
improvement  in  miss  distance  over  that  of  conventional  guidance  exceeds 
two  orders  of  magnitude,  and  gives  virtual  target  contact. 

Table  3.2.  Miss  Distances 


Scenario  2 


Scenario  5 


Pro-nav 

Pursuit 

Reachable  Set 

T  =  .1  s 
s 


3.6.2  Computational  Delay 

The  scope  of  computation  described  in  Sections  3.4  and  3.5  is  clearly 
a  significant  factor  in  implementation.  In  practice  the  length  of  time 
from  measurement  availability  to  completion  of  control  calculation 
ultimately  governs  the  rate  at  which  course  refinements  can  be  made. 
Clearly,  this  implementation  delay  serves  to  degrade  performance  since 
actual  application  of  the  control  effort  comes  when  the  measurements  have 
lost  some  validity.  To  study  the  implications  of  this  effect  for  the 
reachable-set  based  approach,  a  series  of  simulations  were  formulated 


allowing  variations  of  computation  time  lag  from  0  to  ‘iO  msec.  The  update 
interval  as  before  was  fixed  at  100  msec.  Figure  3.1?  shows  these  results 
graphically  for  the  two  scenarios.  The  performance  behaves  in  a  roughly 
parabolic  fashion  as  time  lag  is  increased.  For  Scenario  two,  the  miss 
distance  went  from  .  1  to  3  feet;  for  Scenario  five,  from  .2  to  9  feet. 
This  is  expected  from  an  intuitive  viewpoint,  since  high  crossing  rate 
would  imply  a  faster  obsolescence  of  measurement  data.  Nevertheless, 
accuracy  remained  considerably  better  than  that  of  conventional  guidance 
under  ideal  conditions. 


Figure  3.12. 

Performance  versus  Computation  Delay 
Reachable  Set  Guidance,  Missile  Body  Reference 
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3.6.3  Measurement  Error  Effects 


The  most  significant  requirement  of  the  reachable-set  approach,  aside 
from  computational  hardware,  is  a  set  of  extensive  measurer,  cnts,  and  their 
associated  sensors.  The  necessary  information  includes  relative 
displacement  (target  range  and  displacement  angle)  and  relative  velocity 
(target  speed  and  heading).  The  required  accuracy  of  these  measurements 
dictates  the  sensor  cost  and  complexity.  To  evaluate  the  performance 
sensitivity  with  respect  to  measurement  errors,  a  third  set  of  simulations 
was  conducted. 

1)  Range  Error.  To  each  measurement  of  target  range  a  white 
Gaussian  error  sample  was  added.  The  standard  deviation  of  the  error  was 
specified  as  a  percentage  of  the  actual  instantaneous  range;  the 
distribution  was  truncated  to  give  only  positive  range  measurements. 
Figure  3.13  shows  an  ensemble  mean  over  25  samples  of  performance  versus 
range  error  standard  deviation.  It  can  be  seen  that  moderate  to  severe 
penalties  result  from  random  range  inaccuracies.  For  a  standard  deviation 
of  less  than  75%,  the  performance  is  still  superior  to  conventional 
guidance. 


2)  Speed  Error.  Controlled  inaccuracies  were  incorporated  into 
speed  measurements  in  much  the  same  manner.  Performance  degradation  was 
more  pronounced  than  for  range  errors,  but  the  basic  shape  remains 
unchanged,  as  seen  in  Figure  3.14. 


Miss  Distance  (ft) 


3)  Target  Heading.  The  angular  measurements  necessary  were 
degraded  by  a  white  Gaussian  disturbance  added  to  each  sample.  The 
standard  deviation  of  the  error  was  given  in  absolute  degrees.  As  seen 
from  Figure  3.15,  both  scenarios  were  essentially  equally  affected  by 
errors  in  the  target  heading  angle.  Five  degree  (90  mrad)  standard 
deviation  resulted  in  an  expected  miss  of  25  feet,  a  substantial 
degradation. 

4)  Target  Position.  For  the  implementation  used  in  the 
simulations,  relative  target  position  was  generated  from  a  measurement  of 
azimuth  angle  from  a  reference  axis.  This  second  angle  is  apparently  far 
more  critical  than  the  other.  As  indicated  by  Figure  3.16,  a  smaller 
standard  deviation  of  three  degrees  (55  mrad)  results  in  the  same  mean 
miss  distance  of  25  feet. 

Summarizing  the  observations  from  the  simulations: 

1)  Given  equivalent  noiseless  measurements,  the  reachable  set 
based  guidance  has  potential  for  tremendous  performance  improvements  over 
conventional  guidance.  This  is  due  to  the  reduction  of  terminal 
saturation  effects,  and  the  even  distribution  of  control  effort. 

2)  Actuator  lag  due  to  computation  time  degrades  performance  to 
some  degree,  although  the  effect  may  not  be  serious.  Higher  crossing 
rates  increase  the  severity  of  the  degradation. 


Target  Heading  Error  (degree) 


Figure  3.15 

Performance  versus  Target  Heading  Error 
Reachable  Set  Guidance,  Missile  Body  Reference 


Target  Position  Angle  Error  (degree) 


2 


Figure  3,16 

Performance  versus  Target  Position  Angle  Error 
Reachable  Set  Guidance,  Missile  Body  Reference 


3)  Range  accuracy  is  apparently  not  critical.  Distance  measured 
by  conventional  means  should  be  adequate.  Image  ranging  might  very  well 
be  possible. 

4)  Likewise,  relative  speed  accuracy  is  not  critical. 

5)  On  the  other  hand,  target  heading  requires  measurement 
accuracy  on  the  order  of  a  few  degrees  for  suitable  performance. 

6)  Target  position  azimuth  angle  apparently  required  even  a 
tighter  tolerance,  but  well  within  present  sensor  accuracies. 

The  preliminary  results  described  here  are  merely  intended  to  indicate 
the  sensitivity  of  a  typical  advanced  guidance  technique  to  Independent 
errors.  The  apparent  sensor  accuracy  required  is  stringent.  Of  course, 
an  actual  system  would  be  able  to  incorporate  a  tracking  or  smoothing 
algorithm  in  software,  reducing  a  sensor's  raw  error  substantially,  and 
consequently  enhance  missile  performance. 
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0266 

c 

0267 

TE».(1)  aCOS  (PH  113  AR)  *COS(THFTPR) 

0268 

TEh(2) aSIN(PHIbAK)*COS (THf  T0K) 

0269 

TEP.  (3)*SIN(ThETtR) 

0)270 

c 

w27  l 

c 

ROTATION  BACK 

0272 

c 

0273 

c 

0274 

TEi’.PaTEP  (?) 

«27S 

TEm(2)  aTfc.l  (3) 

0276 

TEPi  (3)  aTEMP 

0277 

CALL  MTRMl(THETM,TfcM,HOLl ) 

H278 

TEnPaHOLO (2) 

0279 

POLO (2) aHOLO (3) 

0)860 

HOLT  (3) a T E P P 

I  "l 


1 1 1 1  T  i '/  r. 


If  . 1  ; 


rnn:F  iGns?  , 


it  ?  M  C  M  l  'If-  I  (f‘i(I  ,  ,1K") 

*?.  a?  C 
/?P3  C 

.264  c  r  «c<  TO  iiPufcrfltM.  IN  ►'ISSIlF  f  Y  ST  £ '  ■ 

L 

»»2fe  6  1 1  *A1  A-  ?(Tr.?*(?)  ,  TF».  (  1)  ) 

l  ?(7  T  f- 1 T  w  s  a  T  /»  k  Ji  (  l  fc  f  f  3  )  #  S.<*«T  (UPT  (T  F  i  •.  f  j)  ,  1 1  •'•*  (1  )  ,  ?  T)  ) 

*2M  C 

tf2t*9  C  M  f»  D  llf-T  V  K  L  T  i  •  A  *  r,  |  ^  Mfl  4  Frl  !■  T  m  1 S  A  L  r  f.L  LN  ft  T  ]  ON.  A^ULL 

2  b  d  C 

/2v  l  CALL  T  I  bT  (k  , » i  ,  ,  f'».  I  i.* ,  T  t  T  D  »  o  I  tj) 

lr  292  C 

f.CCLCC  FImIShKL  AT  T  it  I  ;i  f r  1  ►  T  ****************  t  *********  t  tt  *******  *  * 
t'P94  C 

lr  ?  V  b  C  COnVLPT  THIS  LIFT  V  f-  L  T  i  1  ’«<  f:«r*  ir  I<  L  k  T  IAL  *EFF»-F,'.Ct 

/2!  6  C 

*.?<  /  CALL  MlSt:F  Cl  ,  X  T  f  *»*,  K  AMir  f  »‘l.  Ak  ,  VT  ,PrtT  ,  the  T  A,H«U; »  MU) 

i’2  9H  1 1  i  ( 1 )  *  (  ■ )  S  (  X  r-  ( ii )  )  *  L  T  i>  ( X  N  i.  F )  5 

i  29  9  sSP’CXn  f:0  ) *LPb  ( XN  (6)T 

fci M  UMi)  =  si>,  (X.,  (bT  ) 

H7\'\  LF.L  U)  =~TF  ^(2) 

t;3l  2  L  F  L  ( 2 )  =  T  Lf  (  1  ) 

4  3.  *  LFL(Jl*«, 

/3.  4  CAlL  V.  i-L  (uFL  ,  t  ,/bi.f-T  fL-CT  (t.fil.  ,■  F  L#  31  )  ,  I  F|.,  oT 

S  bK-sAT  A>'2(UC‘I  (Fif,Ln,rM  ,.1),nPLI'  (3)) 

i?'  >i  IF  (A-s, r. 0,-1  ,t  AH  fifi  TP  14P 

V.»>  /  t.  (  1)  sTi  k-.ST 

t'3(  »i  L(21*A.-. 

l-.3i  9  L(,*.)=MG 

^  •  3  J  #  1  4  i/  «  F.  t  II  ri  M 

2311  c.  v(> 


*  *  N  0  fc  r  9  L  F  $  * 
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f  AUfr 


(FTK/1--M  !  £  A5*  ?<!/7C— »ulU  Si,  11*74) 


•• 0  1 


Mf  1  F  T  ,  ( 4  ,  L 

00*2  S  U  ■}  P  »J  li  T  I  NE  H  1S*F  (  TFl.  AG  ,  *T  ,  x"  ,  H  AM?F  ,  1*b  At  ,  vT  ,PnI  ,  ThET  A,  A  VEC  ,  S  JG) 

M£*3  C 

00*4  C  30^711  ROUTINE  T  f.i  cn-.PUTt  lji*.SF*v  ATinr.s  In  hlsSRE  CEMflPfcP 

0005  C  C  0  0  W  D I K  A  T  t  SYSTEM,  GJVFN  INlRTIaL  RtFE’ENCE  COORDINATES 

0M6  c  60*714  MtlO  Fop  KIGhT-PaEl-FD  CPO«uI\ATE  SYSTEM 

0007  C 
K008  C 

ree9  c  1FLAG  =  t.,  GIVE  f-E  TmE  MEASUREMENT,  *  0  kF.IUkn  9  tTK  At  SF  Or  HEP  a  V 

00  1  '4  C 

0011  c  XT  a  TARGET  STATt  VECTOR,  (X,Y,Z)  X  ( V ,  Ph  I ,  T  i.F.  T  A ) 

Ml  2  C 

0013  C  XF  a  M I S  S  I L  F.  STATE  VfcCTr*  hTC 

0014  C 

0015  t  RANGE  *  SEPARATION  riSTAisCt 

0.0  16  C 

0017  C  cEAR  a  AZIMUTH  ANGLE  TO  TAEKET 

00  18  C 

0019  C  VT  t  TARGET  SPEED 

0M20  C 

M21  C  PHI  s  FEARING  ANKLE  OF  TARGET  VFLPCTIY 

01*22  C 

0023  C  |  H£T  A  a  ELEVATION  ANGLE  OF  TARGET  VFLOClIY 

0M24  C 

0025  C  AVF.C  a  VECTOR  FOR  CCMfiOL  EFFORT 

0026  C 

4027  C  SIG  *  A OGLE  FOR  LIFT 

M2*  C 

0029  L  I  M  £  f,  S  1  ft  N  XT(!),Xb(l) 

003*  DIMENSION  b>El(3)  ,  VMVEC1  (3)  ,VKVEC2(J)  »VMVEC3(4)  ,VTVfcC(3) 

0031  DIMENSION  i<  VEC  ( 3 )  ,  HOl.P  1  ( 3 )  ,  HOI  02(3) 

tit‘32  I!rEn$IOn  AVfcC  (3) 

UU33  C 

4034  C  GOING  UR  COMING 

0P35  C 

0046  IF(IFLAG)  m, 9,100 

00  37  C 

0036  C  COMPUTE  RFLATIVE  riSPl  ACEEE^T  IN  I^TERTTaL  rf  PRO  CaU  S 

01*49  C  AND  RANGE 

k,  0  A  \!  c 

et'Al  9  CO  10  K»l,3 

0  *42  10  CEL(A)sXT  (  K  )  »•  X  0  ( tt ) 

004J  RA,*GE*S0kT(MI)T  {  r  E  l  ,CEL,3)) 

00  4  4  C 

0043  c  COrPUTt  VELOCITY  UNIT  VECTORS  IN  INERTIAL  COuRClNATFS 

0.146  C 

4047  VM  VEC)  (1  )  aCPSCXI  (fS))*ru5  (Xi- f6)  ) 

004»J  VPVEL1  (?)  aSln  (X"  (5)  )  *0 05  (XI  (b)) 

0049  VMVEC1  (3)8&In(xi  (tO) 

0(*50  vTvEc  ( l)  scu3  (xt  CO  )  *crs(x  T  (rO  ) 

00  51  vTvEC(2)aSIN(X  I  (b))*C(  Sf*T(fO) 

fc052  VTVfcC(3)*SlN(xT(fc)) 

H«53  C 

0064  C  COMPUTE  RANGE  UMT  VtCTC* 

0*53  C 

0053  CALL  VM!lL(RVtC,l./rtANGE,l  FL.3) 


*  a  i-  f.  ■■■,.?  (F  if /'.--ru.  r  asf  <?  4  i  /  7r>- aijcms  t  ,  u>74) 


ttt-/  L 

l  (  I'*1  C  l  i  '.ST.<Ui:l  "Y"  UMi  vHll-f- 

U^j  L 

i,  if)."  (.  all  V' nl  (M{)Lf'i ,  I  t‘H.  vtr  ,  Vi  Vr.r  1 ,  s)  ,  w  vffh.O 

l ri 1  L*LL  V^lr(VHV^LP,^VpL,nfL[l,i) 

k-.Vf>3  CALL  VMiLtVHVfcC?,  1  ,/Sf'KT  CriiT  (V  V  (-.  C.  2  ,  v  4  F  L  ?,?)),  V  r.  V  fc  C?  ,  0  ) 

i  i  r  i  c 

i?b4  C  CCiSTMCT  "7"  JAMS  VfrCTOP  F  •<(-■*■  X  L«OSS  Y 

M  b  5  V^VtLAU)=Vr,VFl.1(3)*VKV^L?C3)-V!‘v£C1(3)*VMVh‘C?l?J 

HCtS  VKVM-3(P)  sVi'VHll  (?)  *Vt  )-k  *'  ytCl  (1 )  *vf  vFC?  C3J 

^|Y^7  V  f  VK(  3  (3)  =V  4t  cl  (  n  *  V*  VF  L?  (?)  -  VI'  vFL  1  (?)  *  V  v  v  £  L  2  ()  ) 

k'  l  i  f  C 

!■••  »  5  C  ..f.Hl.Tt  FPSIflUH  A  2  I  ‘H  T  H 

**  £  7  *A  L 

d'-<7  l  T  F  •  PsH-JT  (F’VF  C  ,  vf'Vf  1 1 , 3) 

('  vn.  IF  (AttS(Tfct'P)  « G  T  •  1  ( )  ir.jiPsSlG -Jfl  ,  ,UlP) 

l'P7 i  t.-t  AKfsACnjj  (Tfcf’P) 

H  t’  7  4  t 

Kir 7  5  L  1 0  i .  H  U  T  F.  V  F  L  P  I.  IT  Y  A  Z  1 1  U  f  H 

n  S?‘ i  C 

k  /  7  7  X  =  U  t  T  (  V  T  V  F  V. ,  v  i- ■  V  F  C  1 , 3 ) 

l'  V  7  5  Y  s  r  PT  { V  T  V  F  F ,  V  >  V  F  C  ? ,  J  ) 

I'l'/T  as'.  FT  C  v  T  vF  C,  v!  vf  (;.«,?) 

C--  f:  •*  f  F  I  s  A  f  »'  f.p  (  Y  ,  X  ) 

HV-1  IF  (MS(  /).!,  T.l.)  7=5  IP  ?) 

HI?  T  h  [  T  »  s  a  S  1  'j  (  /  ) 

IF  3  V  I  =  X  T  (  u ) 

FI  4  h  l.  T  l '  K  N' 

1-.  iC  3  C 

C  f»E.Titf«N  I  IF  T  V  f.  L  T  t  F 
YV67  C 

,Wi-a  l.'Y  TAIL  V.'ilL  (fiDL  L>  i  ,  blf  (Ml')  ,  Vi  \>r.r?,  3) 

Wl‘h!i  CALL  v-iil  (hOUf?»  CPS  fMU)  »VPVF.C3,3) 

P  l'  L'  •'  (.Alt  V  AIK)  (  AVkC,HlU.Dl  ,HUU?,  3) 

F  t' 1.1  h  £  T  IJK  N 

*'  r*  H  2  t  *■  !  I 


*  *  \(J  F  k  "  (.'  V  * 
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f  AfcF.  :Vi  1 


(PTf«.-MI  t/iSF  ?41  /7C--Ak'P-'S1  *  lu/'l) 


H/El 

HN 

i  L 

U  \'  v  2 

6UHRUU1  !*.•/.  I^ThR  (X  is  K  ,  i:t  L7  A  ,  VT  ,  V  ,  TTr-  ,  Ml)  ,  LP/.i'.nF*  ,  lThK  ,  UK) 

p/c-3 

C 

/Pt  4 

C 

801-532  kOuTlM't  10  S  ,)L  V  F  Fl'P  VELOCITY  OP  I E  *s.  T  Al  ION  FO"  InIEKCEPT  1 

//k  5 

c 

I.JVEn  T  AkGtT  J  OF  C  P  E  A  (  J  00  A  N  0  i  OMML  ^  I  5  £>  T  U  F  S  E  E  1 0 

00/6 

c 

00/7 

c 

0t»  tft25  M)U  TO  USE  ALTEht  ATE  SOUlTIOr,  ALGtL;nAlC  EUM'^llO^ 

OF 

00/8 

c 

AMilfc  vaniarles  fjkst 

00/9 

c 

/>  0 1  / 

c 

X  I  U  a  INITIAL  MJSSILh  POSITION  wPT  TL  TA*r,fcT  (xft,Y/,Z0) 

0011 

c 

(TARGET  INITIAL  VELOCITY  COLll  EAk  TO  X  Axis.) 

0/12 

c 

O013 

L 

K  s  TAKRtT  F'AX  TUHN  kATE 

00)4 

c 

//-15 

c 

OELTA  s  TARGET  TUKN  PLANE  INCLINE 

U  0)6 

c 

00  17 

c 

vT  *  T A P 0 E T  NOMINAL  SPEED 

00  1  6 

c 

0019 

c 

V  a  FISSILE  NDnlNAL  SPEED 

0020 

c 

0021 

c 

TIN  a  INITIAL  VALUE  FOP  GUESS  OF  InTc'PCEPT  Tl^t  (IF  tr, ,  COMPUTES  i 

UZ?2 

c 

0023 

L 

PIT  a  VECTUP  OF  INTERCEPT  INF  (j«K  AT  I ON  ( T  IF Y  ,  F  >■  I  ,  1  PET  A  ) 

002  4 

c 

0025 

c 

EPS  s  h  A  X  In  Lin  KAt-T.fc'  K  F  *  C  it  S M  u  A  k  F I)  TO  TlJLEPAlfc  l  T  ?.PUl  T  F  i>) 

0026 

c 

0027 

c 

f-lTEP  a  MAXI, 00.  -.Ui  PEP  CF  oF-TUN  STEFS  To  Iwy  (I.vPuTlEtO 

IH2  8 

c 

0>0  29 

c 

lTt«  a  MOMuFP  uF  NENTCN  STEPS  Lr...PLF!EU  (OUTFOI  TtO) 

00  3/ 

c 

0/31 

c 

1EP  a  ERPUH  kETOkN  NiICUF.K  a/  OK 

0/32 

c 

*1  SINGULARITY  LPiOITIPn 

//33 

c 

=  -1  FaTL'IKE  TU  LOi  VtFC.Fl)  IN  RJTER 

STH 

0.0  34 

c 

0/35 

c 

0036 

0 I PENSION  XIN(3),F1T(?) 

0/37 

LIi-ENSION  SCf<l  (3)  |SCK2(3) 

0038 

PEAL  K 

/P39 

mATA  PI/3. 1415*27/ 

004/ 

DATA  VNO / 725/0,/ 

00  4  1 

c 

0/4  2 

c 

PANE  INTELLIGENT  INITIAL  GUFSS  (  T  A  11  C^ASt) 

/Z43 

c 

// 44 

V  P  «  V 

0045 

IF  (V,LT#VN()mj  VMsVI.CM 

00  46 

ITfePa/ 

k,Z47 

TFbTIN 

0/48 

IFITIn,GT,«.)  i.O  TU  1 

P/49 

KANGtssnpT (OUT (X  IN , XIU,  ?)  ) 

0  P  b  ft 

TFsRANGF./VM 

//51 

c 

0/52 

c 

NEwTUN  STEP  LOUP 

0/53 

c 

0/54 

1  >  s  (V«*TF)**?-Cxlf  (D-VT/c  *SIn  (n*TF)  )  **2 

tf/55 

*  “  C * I F  C2)+VT/KaC(iS  (CELT <)  *  (COS  (K*TF  1  -1  ,)  )  **P 

0/56 

fc  -  (XIN (3)+vT/K»51N (TtL  T  A) * (OUS (F  *TF) -1 ,)) **P 

63 

I’ti.r  ,  *;■  I»'TF*  (F  TN«:--Kfc|.  c*sf  2  4  1  7  7  C- - :  r.  i.s  f  ,  1<70 


l  ► 

p  P,  ts  p 
p,FbB 
t'f  t>7 
tiH'l 
0(-f>2 
0  ?  i'  3 
O.’o4 

0  i"  f>  b 

f  t  c  7 

00  f’  Q 
v‘  i'  7  i", 

02  71 
0*72 
t'k73 
0  t  i  A 

k‘1'7  9 
00  7b 
00  7  7 
0k  78 
H70 
*-i0(7 
7  0  l'  1 
1-0  r  2 
Vi  t  3 

P  »  f.  4 
0  0*5 
00o6 
7,787 

fc.0  »8 
0  0  o  0 
t  2  .7  / 
lr)0  N  1 

v  "  y  2 
0  0  !«  3 


tn  =  ?,*V'*V:*fF-?,*(XlN(l)-vT/0*fJIK(K*TM)*(-VT*ru?  *TH) 

*  - ,] ,  *  (  V  I ( 2 )  ■»  v  T  / ,.  *  U  is  C  i  FI  T*wrcCb(**K')-l.n 

*  *  (>VT*rijS  (Hhl.1  A)  »SIi.  ('  tim 

>.  -2.*f*T»l(i)+V7/K*SlK(  ;Ll4l*fcro(8*1Fj-l.))*f"Vl*Mr4(ljFLTA) 

7.  *  S  I  ■  (  K  *  1  F  J  ) 

IF  (  a i •  s  ( v i ' )  ,(.1 1  j  a t  ~ .‘j )  r.Li  1  r  P-) 

I  F  (  S  U  V  T  f  *•  0  s  (  Y  )  )  .  I.  E  ,  t  F'  .c  *  2  r  *  T  F  )  l.,n  Tu  F  9 
rFsTF-Y  /  (!' 

IF  lIF  ,1  I  IFs.M 

x  T  K  <  a  I  1  fc  <♦ 1 

IF  l  Ilf  *  .fit  ,”ITt8)  M*  U.  1 1-  9 
l«  0  T  fj  1 

(. 

C  C  k  M.Tol''< 

C 

F  S  I  E  (•  s  / 

*«C.a  t-Kju  {3)-VT/K*SIf4(IFLT  Al  *  (LC3  -1  .  )  )  /  7V"*TF  1 

]F  (AmS(  AKK)  ,GT,1  .)  A  MG  *S  I  f  i.  (I  .  ,  A  *G  ) 
ihctasasika.u;) 

(-rf  JKI  (g)-Vl/K*Cn.Sirtn  A)  *  (CuS  (0*TF  )-1  ,n 

K  i;  f  I.  s  (  .  X  I  N'  (  1  )  +  V  1  /  k  *  5 1  ‘1  (  K  *Y  r  )  ) 

FhisaI  At 
«S  HlTfiJsTF 
FIT (0] s  p  h I 
dl ( jJslnFTA 
K  t  T  I  ’  K  N 
C 

C  i  I  ■  S >- i L  /» Ft  I  7  y  uf  ri 

C 

y*  ifksi 

I.F1UKW 

L 

C  F  A  i  I. li ^ L  Til  CUf.  YFKpF  M  Tl.KK 

L 

1  0  '■>  IK*.  =  -1 
m  F  1  1 1  h'  l\ 

I.  M 1 


I 


*  * 


\f.)  F.rPM-b* 


a>  n  a  oi  &  onu 


F-'ACfc  il-.  1 


(Fir  /I  —  f.  1 1  CASE  JU1 /7C--AURURT  ,  ly7<0 


f  T  v  ^  ,  L 

<?  M'ORUUT  TnE  M1KaK(AKG,VIn,  Vimn 

•J  c 

4  L  *7,711  KOUTI^t  TO  T  K  A  i*  S  f  T  h  ►'  In  RE  1  A  ]  Ii-n  aI  SENS*7  *  3-SFACE 

5  C  VECTOR 

b  C 

7  C  ANG  b  ROTATIONAL  ANGLE 

8  C  " 

$  L  Xlh  s  3-SPACt  VECTOR 

*  C 

I  C  X OUT  8  ROTATED  vector 
c 

LlMEuSlOi.  X I N  ( Ck 5  ,XUUT  (d) 

*CUT(3)«XlN(a) 

XOuT(l)aCPS(AM.;)*XIK  (  1  )  ♦  i>  I K  f  AF  G  )  *  X  U  (  ?  ) 

XOUT  (5>)  s-SIN  (AirP)  * X  I  r  (1)  +CU5  (  AnP)  *XTn  U) 

RETURN. 

END 


NO  ERRCRS* 


i  t  1;  t  1 


It  P  I  r.«SF 


M; f  k i  i •  T  I ;  •  F  <1  J  (  Ai.n,  X  I*  ,  yi.M  n 

’  c 

I  C  a  1  V  K  •'  ii  t  T  A  t-  S  F  ( i  fc  i .  A  T  1 1 1  ► 

I  L 

I  1. 1 '.FI;  SI  (If.  /I,-(l)  ,y.]l  Till 

»  X  C  li  T  (  i  )  s  X  I N  f .»  ) 

>PuT(l  js(f:i( 

XnuT  tP|  f  X  IK  Cl)*crb  f  AMil  *XIN  l») 

i  >t  F  T  l .  k  v 

r.  \i; 


M  uS  i 


F  0  t  F  K  l  a  s  * 


